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ABSTRACT 

We investigate how environmental effects by gas stripping alter the growth 
of a super massive black hole (SMBH) and its host galaxy evolution, by means 
of ID hydrodynamical simulations that include both mechanical and radiative 
AGN feedback effects. By changing the truncation radius of the gas distribution 
(-Ri), beyond which gas stripping is assumed to be effective, we simulate possible 
environments for satellite and central galaxies in galaxy clusters and groups. 
The continuous escape of gas outside the truncation radius strongly suppresses 
star formation, while the growth of the SMBH is less affected by gas stripping 
because the SMBH accretion is primarily ruled by the density of the central 
region. As we allow for increasing environmental effects - the truncation radius 
decreasing from about 410 to 50 kpc - we find that the final SMBH mass declines 
from about 10^ to 8 x 10® M©, but the outflowing mass is roughly constant at 
about 2 X 10^° Mq. There are larger changes in the mass of stars formed, which 
declines from about 2 x 10^" to 2 x 10^ Mq, and the final thermal X-ray gas, 
which declines from about 10^ to 5 x 10® M©, with increasing environmental 
stripping. Most dramatic is the decline in the total time that the objects would 
be seen as quasars, which declines from 52 Myr (for Rt = 377 kpc) to 7.9 Myr 
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(for Rt = 51 kpc). The typical case might be interpreted as a red and dead 
galaxy having episodic cooling flows followed by AGN feedback effects resulting 
in temporary transitions of the overall galaxy color from red to green or to blue, 
with (cluster) central galaxies spending a much larger fraction of their time in the 
elevated state than do satellite galaxies. Our results imply that various scaling 
relations for elliptical galaxies, in particular, the mass ratio between the SMBH 
and its host galaxy, can have dispersions due to environmental effects such as 
gas stripping. In addition, the simulations also suggest that the increase in AGN 
fraction in high-redshift galaxy clusters might be related to environmental effects 
which shut down the SMBH mass accretion in satellite galaxies and reduce their 
AGN activity. 

Subject headings: galaxies: active — galaxies: elliptical and lenticular — galaxies: 
ISM — galaxies: nuclei 



1. Introduction 



The role of the environment in galaxy evolution has been suggested in various forms 
which strip out g as fro m a galaxy and goes all the way back to the early suggestion by 



Spitzer fc Baadd (119511 ). For example, the ram pressure of the intracluster medium is a 



ters and to stop the supplv of cold gas for star formation ( Gunn & Gott 


1972; Larson et al. 


1980; 


Takeda et al. 


1984; 


Gaetz et al. 1987; E 


5egelman & Fabian 199ol: 


Abadi et al. 1999; 


Domainko et al. 


2006: Tonnesen & Brvan 2008 


; KaDferer et al. 2009). Ot 


ler possible pro- 


cesses include thermal evaporation and viscous s 


tripping (Cowie & Songailal 


I977I Livio et al. 


1980; 


Nulsen 198 


2; 


N^epveu 


1985; Valluri & Joe 


1990: Roediger & Briiggen 


2008|). The com- 



bined g as loss by these different types o : 



groups (iStevens et al.lll999l:lQuilis et al. 



20081 : iMcCarthv et al 



2008 



destruction effects is expected in galaxy clusters or 



Smith et al. 



2OOOI: Toniazzo fc Schindlerll200ll:lKawata fc Mulchaey 



2OO9I ). Tidal stripping can also p lay an impor 



tant role in changing the gas and stell ar mass of cluster or group galaxies (IMerrittl Il983 
Moore et al.]ll999l : lo'Ercole et al.lboooh . 



Recent multi-wavelength observations have proved the loss of gases in cluster or group 
elliptical galaxies. In some elliptical galaxies, diffuse X-ray emitting gases show a long 



Kim et al. 2008a 



tail structure which can be explained by ram pressure stripping (e.g 
Randall et al.l 120081 ). Infrared observati ons also revealed d ust emissions that can be from 
gas stripped from elliptical galaxies (e.g. IWhite et al.lll99l[ ). 
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A central galaxy which hosts satellite galaxies in groups and clusters exhibits different 
features compared to satellites which can be affected by the various stripping processes. 
Because the central galaxy sits near the bottom of a deep gravitatio nal potential well, its hot 



Sun et al. 20051), and it is not surprising 



gas halo is much larger than those of satellites (e.g. 
that it shows signatures of cooling flows in some cases (see lFabianlll994j : iReiprich et al.ll2004j . 
for a review). 

Quite obviously, the differences in the stellar populations between satellites and central 
galaxies also have been studied, being relevant to our understanding of the environmental 
effects . For example, red satellite galaxies are redder than central galaxies of the same stellar 
mass (Ivan den Bosch et al.ll2008l ). Many central galaxies in clusters display either recent star 
formation or ongoing sta.r forn i ation which may be related to the cooling flow from their hot 



gas halos fICardiel et al.l Il998l : iRaffertv et al.l 120061 : iBildfell et al.l 120081 : lO'Dea et al.l 12008 



Pipino et al. 20091 ) and which would be of reduced significance for satellite galaxies. Satellite 
galaxies in galaxy clusters seem to gradually lose of their gas by environmental effects and 
to truncate star f ormation, cornparin g their current and past star formation with those of 
field galaxies (e.g. iBalogh et al.lll999l ). 



Differences in the properties of a central super massive black hole (SMBH) and its 
activity are expected between central and satellite galaxies, when considering the strong 
correlation between spheroidal galaxies and their SMBHs. The ratio of the SMBH mass over 



in local galaxies ( 


^ormendv & Richstone 


I995I: 


Maeorrian et al. 


19981: 


Ferrarese & Merritt 


2000 




Gebhardt et al. 


2000: 


Ferraresd 


2OO2I: 


Yu & Tremaind 


2OO2I: 


Marconi & Hunt 


I2OO3: 


Harine & Rix 


200^ 


1). If the growth of stellar mass and SMBH mass are coupled to the same 



kind of environmental effects, it is natural to conclude that the properties of SMBHs in satel- 
lite galaxies must be somewhat different from those of central galaxies. As we emphasized in 
our previous papers (Ciotti et al. 2009a, hereafter Paper I; Shin et al. 2010, hereafter Paper 
II) the self-regulated growth of SMBHs is profoundly dependent on how the accretion energy 
is converted to heat the ambient interstellar medium and how frequently and quickly the 
feedback process is ignited in the right place: we called these two key issues the problem of 
energy conversion and the timing problem. One would expect that the environmental effects 
might induce differences in the frequency of AGN feedback and the growth of the SMBHs 
by varying the conditions of the self-regulation process. 



primary galaxy accreting satellites in group and cluster environment corresponds to the brightest 
cluster galaxy (BCG) in this paper. For isolated galaxie s, the term central galaxy is correspondent to a 
primary galaxy which hosts satellites (e.g. lAnn et al.1120081 ). In this paper, the key feature to define satellites 
is that they experience gas stripping which is strong enough to affect their evolution. 
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In this paper, we tackle the issue of environmental effects on the coevolution of the 
SMBH and its host galaxy by simulating the evolution of a galaxy in hydro dynamical mod- 
els with a simplified setup of different environments. Although there is no current systematic 
investigation of the SMBH mass and AGN activity for separating satellite and central galax- 
ies, the theoretical prediction from our simulations can be used to constrain the hypothesis 
of the coevolving SMBH and its host galaxy and the connection to environmental effects. 

This paper is organized as follows. In Section [21 we describe the models and the simu- 
lations. The results are presented in Section |21 Discussion and conclusions follow in Section 

SI 



2. Simulations 



We adopt as a basic model the BS"^ model that is given in ICiotti et al.l ( l2010l . hereafter 
Paper HI). This model incorporates both radiative and mechanical feedback effects from the 
central SMBH in ID hydrodynamic simulations. The basic input physics is fully described 
in Paper I and II. Here, we summarize the main ingredients of simulations and explain how 
to implement different environments. 

We setup an initial galaxy following observed properties of local elliptical galaxies. First, 
the stellar mass is distributed to be consistent with Faber- Jackson relation and the Funda- 
mental Plane, assumi ng an identical amo unt of stellar and dark matter mass within the 
half- mass radius (e.g., iTortora et al.l 120091 ): all the relevant dynamical and structural prop- 
erties of the models are given in Paper HI. Second, the initial mass of the central SMBH 
Mbh is 0.001 of the initial stellar mass M*. In all simulations, we set = 2.9 x 10^^ M©, 
the effective radius Re is 6.9 kpc, and the central aperture velocity dispersion is 260 km/s 
(see Appendix for different initial velocity dispersions and their effects). Therefore, in our 
simulations, we are tracing the late evolution of ellipticals after the initial epoch of formation, 
i.e., of major growth of both the stellar component and of the central SMBH. 

In addition to passive stellar feedback effects, i.e. mass losses and type la supernovae, 
the simulations include both radiative and mechanical feedback due to the SMBH as well 
as secondary stellar feedback effects from recurrent star formation. As discussed in Paper I, 
II, and HI, the observational constraints such as the mass ratio between the SMBH and its 
host galaxy, the X-ray luminosity from hot diffuse gas, the quasar lifetime, and the recent 
star formation fraction, are better matched by simulations which include both feedback 
mod es than simulati ons restricted to purely radiative or purely mechanical feedback models 
(see iBegelmanl l2004j . for a review of various forms of feedback). Radiative feedback affects 
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the ambient medium around the SMBH via radiative heating mediated by photoionization, 
Compton heating, and radiation pressure. These can be effective forms o f AGN feedback over 
a lar ge spatial range depending on a frequency r ange of AGN radiation ( ISazonov et al.ll2004j . 



2005 



Fabian et al. 



2006t ICiotti fc Ostrikerll2007( ). Meanwhile, the typical form of mechanical 



feedback is represented by the interaction of nuclear winds and jets from accreting SMBHs 
which ultimatel y heat up the surrounding interstellar medium with maximum effectiveness 



in r <C 0.1 kpc (jTabor fc Binnev 



20051 : ISpringel et al. 



2005a 



1993 



Konigl 



2006 



Johansson et al. 



Binnev fc Taborlll995l:lBegelmanll2004l : IVeilleux et al. 



20091). 



In our models, we control the strength of radiative feedback with the parameter eo which 
is related to the radiative efficiency esM that determines the bolometric accretion luminosity 
from the SMBH 

^BH = CemMbhC^ (1) 



by 



Atti 



l + Am 

Equation [2] is a phenomenological implement ation of the main fea tures of the ADAF (Ad- 
vection Dominated Accretion Flow) accretion ( iNarayan et al.lll996l ) with the rescaled SMBH 
accretion rate with respect to the Eddington accretion rate rfi = M^u/ M^dd- The parameter 
A is 100 as in Paper I and II. We use eo = 0.1, which is favored b y observational constraints 
for luminous quasar phase, i.e. when m is high (e.g. ISoltanlll982[ ) §. 



The parameter handles the mechanical efficiency of the AGN wind which is coupled 
to the bolometric accretion l uminosity of AGN and affects the mass loss rate of gas by the 
wind (IKurosawa et al.ll2009al ). After the introduction of the normalized accretion luminosity 
with respect to the Eddington luminosity LEdd, 



I = 



^BH 



Atti^ 



Lsdd 1 + Am ' 
the mechanical efficiency of the AGN wind e„ is 

3e!^ / 



(3) 



4 1 + 0.25/ 



(4) 



In this paper, we fix = 3 x 10 ^ which is one of the explored models and suggested as the 
probable value in Paper HI. The peak mechanical feedback efficiency governs how much 



^In Paper III, we also present models with eq ~ 0.2. 



- 6 - 



kinetic energy, momentum, and mass h as to be deposited to the b road-line region, with the 



mass loss coefficient of the AGN wind (iKurosawa fc Proga]l2009bl Jd) 



4 1 + 0.25/- 



where 77^ = 1800e^ which is consistent with observational constraint s on the velocity of th e 



broad-line regions and the ionized absorption outflows in quasars (e.g. iKrongold et al.ll2007l ). 



The important feature of our simulations is the boundary condition that mimics different 
environments having different strengths of gas stripping. In practice, this is accomplished by 
changing the position of the last grid point Rt (hereafter called "truncation radius") where 
outflow boundary conditions are imposed. We assume that gas escapes beyond the truncation 
radius as various destruction processes caused by the environment (pri ncipally, ram-pressure 



stripping) reduce the size of hot gas halo, but do not totally destroy it ( iBrighenti fc Mathews 



I999I : ISun et al.ll2005l ). There is no infall of gas from outside these radii. 



The e fficiency of g as stripping depends on both galaxy and intracluster medium prop- 



erties (see 



Hester! l2006l for discussions). The gravity from mass distribution in the galaxy 
competes with ram pressure which depends on the density of intracluster medium and relative 
velocity of galaxies with respect to intracluster medium. Therefore, the stripping strength 
can vary significantly depending on the complex combinations of galaxy models and the flows 
of intracluster medium toward galaxies, i.e. the orbits of galaxies and the properties of the 
intracluster medium. 

Even though this setup of different environments is too simple and the symmetric strip- 
ping never happens in real conditions, it will help us to understand how environmental effects 
might affect galaxy evolution with AGN feedback. Simply, small values of Rt correspond to 
the satellite-like environment where gas stripping can be effective, while large values of Rt 
represent the environment of central galaxies which can keep a large amount of diffuse hot 
gas, practically without any substantial loss of gas (but without the confining effect of some 
external pressure, as for example, the intracluster medium of the center of a galaxy cluster). 

We test ten different truncation radii for the same initial stellar mass and SMBH mass. 
As given in Table [H Rt ranges from about 51 kpc to 413 kpc. This range corresponds to 
about 7 to 60 Re- Initially, all galaxies have the same density and temperature distribution 
of gas within the truncation radius. Therefore, when density and temperature distribution 
of hot gas is same in all cases at a certain time, the total X-ray luminosity emitted from 
the hot gas is expected to be lowest for the smallest truncation radius, i.e. Run 1. The 
spatial resolution of all runs is exactly the same despite the differences of the truncation 
radii, having different numbers of radial grids. The inner most grid point is placed at 2.5 pc 
from the center in all runs. 
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The simulation follows the evolution of a galaxy, which begins at the cosmic age 2 Gyr 
and stops at 14 Gyr, by solving the hydrodynamic equations of the gas, being supplemented 
by prescriptions for star formation and stellar feedback. Because the simulations begin at 
2 Gyr, i.e. a redshift of 2; ~ 3.2 for the LCDM cosmology with = 0.3, Qx = 0.7, and 
Hq = 70 km/s/M pc, it is assumed that a bulk of stellar mass is already established at the 
beginning epoch ( lRenzinill2006l : Ivan Dokkum et al.ll2010[ ). 

As explained in our previous papers, the simulations are not meant to reproduce the 
cosmological context. Instead, we explore the physical link between the local small-scale 
physics around SMBHs and the global scale of a single galaxy during late evolutionary 
phases. In general, it is already difficult to achieve this dynamic range with the intricate 
prescriptions of AGN feedback, so that c osmological gas infall o r galaxy merger/ac cretion 
are not implemented in our simulations ( iDi Matteo et al.l l2008l : iLetawe et al.ll2008l ). This 
limitation must be considered when interpreting our simulation results. However, the recent 
observations suggest that the final stage of fueling onto the SMBH is controlled primarily 
by the self-regulation process instead of large-scale episodic effects such a s galaxy mergers 



(IGrogin et al. 



20091 : 
only. 



2005 



KoUmeier et al. 2006: Li et al. 2008: Gabor et al. 2009: Silverman et al. 



Reichard et al.l 120091 ) . Moreover, in this paper we focus on effects by gas stripping 



3. Results 



3.1. Evolution of models 



In all computed models, the growth of the new stars and the SMBH is strongly regulated 
by how frequently active phases of SMBHs occur. At early times, all models are characterized 
by the accumulation of recycled gas and its cooling which finally results in the formation 
of new stars and almost simultaneous SMBH growth, with a strongly intermittent activity 
as shown in Figure |T] The coincident activity of star formation and SMBH accretion is 
consistent with blue colors of quasar- hosting elliptical galaxies at low and high redshi fts and 



the recent star formation in AGN hosts (e.g. iJahnke et al.ll2004j : ISanchez et al.l 120041 ) 



The repetition of star formation episodes and the active SMBH phases is basically caused 
by self-regulation of AGN feedback, which heats up the cooling flow for star formation and 
suppresses fueling to the SMBH responding the cooling rate. The time intervals between the 
peaks of star formation rate M^, and the central SMBH growth rate Mbh become progressively 
longer and longer, even though the precise timing of the repetition varies in each run. This 
delay is the result of the decline of M^, with cosmic time, and of the competition between 
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heating and cooling which finally resurrect star formation and fueling onto SMBHs. The early 
frequent activities of forming stars and feeding SMBHs result in the substantial increase in 
stellar mass and SMBH mass within the first two Gyrs as recycled gas from stars is used to 
form new stars in the central dense regions. 

The dominance of SMBH growth over star formation changes over time, depending on 
how quickly feedback from the accreting mass affects the ambient interstellar medium (ISM) 
around the SMBH. For example, we investigate the evolution of Mbh and in detail around 
3.23 and 13.31 Gyr in Run 6. As shown in Figure |21 the peaks of Mbh in the early accretion 
phase exceed the Eddington accretion rate for some short period of time. This high accretion 
onto the SMBH finally produces a strong feedback effect which suppresses the star formation 
rate. The late SMBH accretion rate, however, is much lower than the Eddington accretion 
rate even in the peak phase, while M^, rises to lO^M0/yr for a moment. The ratio of 
to Mbh in this late accret ion is consistent with the observed ranges of low-redshift quasars 



as shown in Figure |2] (e.g. IShi et al.ll2009l ). In the late accretion, the ratio of M* over Mbh 
ranges from about 10 to 300 for most time when Mbh > IMQ/yr. This range is consistent 
with the observed range even though the observed range might be limited by samples of 
quasars. Further discussion of this feature will be given in the next section. 

As shown in Figure [3l the mass of the hot. X-ray emitting gas in a galaxy does not change 
as much as the increase of the newly formed stellar mass or SMBH mass. We measure the 
total amount of the hot ISM within lORe in all models except for Run 1 (where Rt is smaller 
than lORe). The mass of outflowed gas is also measured at lORg. In the case of Run 1, these 
two quantities are instead measured at Rt. In general, stronger evolution is found in models 
with a small truncation radius. This has important consequences: for example, in Run 1 
and 2, the decrease in gas mass is so substantial that intensive star formation and SMBH 
accretion does not resume, even though heating by stellar and AGN feedback effects is not 
comparable to that in Run 9 and 10. 

Importantly, the total amount of outflowing gas depends only weakly on the truncation 
radius (see Figure |3]). Despite the small Rt, in Run 1 and 2 the total amount of the gas 
blown out is not higher than those of other runs because they have less frequent intensive star 
formation and SMBH accretion events which dump out energy and momentum to produce 
outflows. We examine this feature further in ^3.31 



Figure m represents the change of the mass budget as a function of Rt at 4, 8, and 14 
Gyr. The difference of the stellar mass appears at 4 Gyr, while other masses deviate less 
from each other. But as galaxies evolve, galaxies with smaller Rt lack more X-ray emitting 
gas and stellar mass. It is naively true that the smaller galaxies have a smaller amount of 
gas and stellar mass because of their smaller Rt. But this trend is not very effective for the 
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growth of SMBHs as we see in Figure HI resulting in dispersions of the ratio between stellar 
mass and SMBH mass. Mbh/^* in Run 1 is about two times smaller than that in Run 10. 



3.2. AGN activity and feedback 



The frequency of AGN activity varies systematically with Rt as shown in Figure [TJ 
Smaller Rt results less frequently in strong SMBH accretion events. Particularly, this change 
is significant for Runs 1 to 3. For example, Run 1 stops any further intensive accretion after 4 
Gyr as shown in Figure O Only early peaks of the SMBH optical luminosity are higher than 
the 10% of the Eddington luminosity. However, Run 6 continues to shows AGN activity until 
14 Gyr, while the late optical emission from the central SMBH is much lower than the 10% of 
the Eddington luminosity as in Run 1 (see Figure [5]). In other runs such as Run 6 except Runs 
1, 2, and 3, this late continuous AGN activity smooths out the dependence of AGN activity 
on Rt, and instead shows highly time- dependent stochastic effects. Remarkably, as we see in 
our simulations, recent observations show a general decrease in the rat io between the SMBH 
accretion ra.te an d the Eddington rate as the redshift decreases (e.g. iKoUmeier et al.l 12006 



Ballo et al.l 120071 ). In sum, AGN activity in satellite galaxies is expected to decrease with 
decreasing redshift much more than that of central galaxies. 

We note that the late SMBH accretion can reach the limit of the Eddington accretion 
rate in models with only mechanical feedback effects (see Paper I and II). But in the current 
models with both radiative and mechanical feedback, the radiative pressure contribution due 
to du st opacity makes the la te peaks of the accretion rate much lower than the Eddington 
rate (jCiotti fc Ostrikerl 120071 ) i n agreement with re cent measurements of the accretion rate 
and hydrogen column density (IFabian et al.l 120091 ). By 3.5 Gyr mass loss from the initial 
bulge accumulates an amount of cold gas which can be strongly affected by radiative feedback 
effects. Therefore, after about 3.5 Gyr all models do not have super-Eddington accretion 
phases as shown in Figure [TJ 

Strong AGN activity is generally correlated with the increase in X-ray luminosity from 
the hot ISM. Despite the differences in Rt and the consequential difference in AGN activity, 
the response of the X-ray luminosity to the AGN activity is similar in Runs 2 and 6 as 
shown in Figure [H As the SMBH accretion rate increases, the oscillatory response of the hot 
ISM to AGN feedback begins to produce an oscillation of the X-ray luminosity. When the 
optical luminosity from the SMBH a ccretion temporally d rops because of extinction, the X- 
ray luminosity quickly increases (see iPellegrini et al.ll2009l . for a discussion on the evolution 
of the nuclear X-ray emission). 
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The effects from AGN feedback can also be found in the radial distribution of gas 
temperature and density. Figure [6] shows the two cases of Run 1 and 6 at 3 and 14 Gyr. As 
we found in Figured both models are in an active phase at 3 Gyr. Basically, the active phase 
in both models causes the increase in the central temperature. In addition, the consecutive 
active phases leave the wiggle structure in both temperature and density distribution as the 
AGN feedback effects propagate to the outer regions on a sound-crossing time scale. At 14 
Gyr, Run 1 is in the quite stationary, hot accretion phase, while Run 6 experienced a recent 
violent nuclear burst within 1 Gyr (see Figure E]). Considering uncertainties of deprojection 
and radial binning effects in observation, the radial temperature and densit y distribution of 



Run 1 is within the observed range of local quiescent elliptical galaxies (see iFukazawa et al. 
20061 : [Humphrey et al.l l2006l . for examples of observation), although its X-ray luminosity 



from the hot gas is lower than the observed range of luminosities for the same stellar mass. 
For example, the central density of the hot ISM is around 10~^^ to 10~^^(gcm~^) in all runs. 
But the recent AGN activity makes Run 6 have a more disturbed radial structure than Run 
1. Moreover, we find a significant difference in the central density between the two runs at 
14 Gyr. 

Finally, the episodic and net quasar lifetime is measured in all runs, defining the quasar 
phase with the optical lumino sity from the SMBH in i?-band Mb < —23 mag (in Vega 



magnitu de) (e.g., iMartini |2004[ ) by using the typical spectral energy distribution of quasar in 
S-band (lElvis et al.lll994l ). This magnitude cut correspon ds to about Lrh 'y 5 x lO^'^erg s~^ 



with the median spectral energy distribution of quasars ( lElvis et al.lll994[ ). We define the 
quasar phase with this Lbh cut and effective bolometric luminosity from simulations. The 
effective bolometric luminosity represents bolometric luminosity from AGN with extinction 
by ISM as explained in Paper I. The episodic lifetime of quasar phase is an interval measured 
from the increase of the SMBH optical luminosity above Mb < —23 to the drop below this 
level, while the net quasar lifetime is the sum of these episodic lifetimes. Therefore, each 
peak of Mbh has its episodic lifetime when it is above Mb < —23. In Figure [TJ we present 
results for Runs 1, 4, 6, and 9. The episodic lifetime of each quasar phase varies significantly 
among the peaks of the SMBH accretion. In particular, an extremely short episodic quasar 
phase is found between long episodic bursts after about 3.5 Gyr. The episodic quasar hfetime 
ranges from about 0.1 Myr up to about 10 Myr. Even though the simulations with the large 
Rt do not experience late luminous quasar phases within the last 2 Gyr (see Figured]), the 
net hfetime of quasar phase generally increases as the Rt increases. But, because the net 
lifetime in models do not include any possible active phase before 2 Gyr, i.e. the starting 
epoch of our models, a direct comparison of this net lifetime to observational constraints is 
not straightforward. 
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3.3. Outflowing gas 

Even though the boundary condition of our simulation might not be reahstic in detail, 
it is useful to know the properties of outflowing gas. Probably, the gas blown out of the 
truncation radius of satellite galaxies will be mixed into the intergalactic medium although 
the gas will change its temperature and density after mixing with the intergalactic medium. 
As explained in §3.11 the outflowing mass is measured at Rt for Run 1, while in other runs 
the mass flowed out at 10i?e, i-e. 69 kpc, is measured as the outflowing mass. As shown in 
Figure [Sj the total amount of outflowing gas does not show a significant difference between 
the case of satellites and central galaxies, while the growth of stellar mass is quite low in 
Runs 1 and 2. Because models with a large Rt do not lose much gas, the outflowing mass 
measured at lORe represents the pollution of the central gas to outer regions where satellite 
galaxies can also contribute gas by stripping processes. Therefore, galaxies with a small Rt 
are more efficient for enriching the intergalactic medium per galaxy. 

Figure [H] illustrates how the high radial velocity and density of gas in models affects the 
total mass of outflowed gas. The galaxy with the smallest Rt (Run 1) has expelled high- 
velocity dense gas continuously. Therefore, despite its smaller outflowing surface area at 
Rt = 51.4 kpc than other runs, its total mass of outflowed gas is substantial and comparable 
to the mass measured at lO-Rg in other runs. For example, the outflowing velocity for 
Rt = 51.4 kpc (Run 1) is close to 400 km/s at 4 Gyr. At 14 Gyr, simulations from Run 3 
to 10 produce the low- velocity outflow which is slower than 100 km/s. But the difference in 
density becomes small at late time such as at 14 Gyr, as we also found in Figure [61 



3.4. Color of galaxies 



Color is one of the well observed properties of elliptical galaxies which are dominated by 
old stellar populations in the local universe. The simple old stellar population of elliptical 
galaxies implies that their star formation histories hav e been significantly qui e t although a 



small amount of recent star formation is possible (e.g. iLucero fc YoungI 120071 : lYoung et al. 



20091 ) . But environmental effects suc h as ram-pressure st r ipping are expected to make the 



color of satellite galaxies redder (e.g. [Martinez et al.ll2008l : IWeinmann et al.ll2009l ). 



In order to investigate the spectral properties of our galaxy models, we construct a color- 
magnitude diagram of model galaxies by synthesizing spectra based on the star formation 
history of the simulations. Assuming solar metallicity or half sola r metallicity for all stella r 
populations, the spectra are synthesized by using the BC03 model (IBruzual fc Charlotll2003l ) . 
and the corresponding magnitudes are estimated from the spectra with the SDSS filter bands 
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( iFukugita et al.lll996l ). Because we do not consider the complicated metallicity distribution 
of stellar populations and dust extinction, the derived color cannot be directly compared with 
the observed colors of elliptical galaxies in various environments. Moreover, because colors 
of old stellar populations fo und in elliptic als are strongly affected by age effects with age 
- metallicity degeneracy (see iRenzinil |2006| . for a review), a direct comparison between our 
simulations with observed ellipticals need the careful consideration of age and metallicity 
distributions. Here, we just focus on how different star formation histories in different 
environments affect the color of galaxies in all runs. For the initial stellar mass at the 
beginning of simulations, we assumed that the same star formation history before 2 Gyr is 
described by 

M. oc (t/r^) e(-*/"\ (6) 

where t is a cosmic time and r 



1 Gyr (e.g.. iNagamine et al.ll2006[ ). Therefore, the peak of 



the initial starburst is at t 



T. 



In Figure [H we present the distribution of model galaxies in the SDSS (m — r) rest-frame 
color and r-band absolute magnitude in AB magnitude at 8.5 and 14 Gyr. All galaxies are 
brighter and bluer at 8.5 Gyr than at 14 Gyr because most of the stellar mass is dominated 
by the initial old stellar population which changes its color by passive evolution. But there 
is a significant variation of the color for simulations with different truncation radii. For 
instance. Runs 8 and 9 at 8.5 Gyr show bluer color than other models, while becoming green 
galaxies temporally. At 14 Gyr, the color of the galaxy in Run 6 deviates strongly from those 
of other runs toward green color because of the recent star formation (see Figured]). This 
temporary transition to green color with the increase in stellar mass can be a reason why 
moderate-luminosity AGN, i.e . a low-accretion AGN, are often fouri d with a green color in 
the color-magnitude diagram (IHasingeii l2008l : ISchawinski et al.l 120091 ) . 



In addition to the stellar flux, the color of galaxies can also be affected by AGN radiation 
if the latter is comparable or more luminous than the radiation from the host galaxy. To 
study this problem, we add the typical spectral energy distribution of type-I luminous quasars 
(IVanden Berk et al.ll200ll ) to the synthesized stellar spectra of model galaxies. The quasar 
spectrum is scaled to the extinction-corrected optical luminosity from simulations, which 
correspond to the observed luminosity of quasars, before it is added. Because the quasar 
radiation is substantially bluer than the stellar radiation of quasar hosts, the addition of 
quasar radiation makes the global color of simulated galaxies bluer, as we show in Figure 
[TUl Particularly, Run 9 is strongly affected by its quasar radiation. The almost coincident 
activity of star formation with SMBH accretion as shown in Figures [T] and [21 results in 
the shift of galaxy colors to the blue by radiation from young stars as well as that from 



quasar phase. Importantly, as observed in A GN host galaxies (e.g lCid Fernandes et al.ll2001 



Ammons et al.ll2009l : [Silverman et al.l 120091 ). the existence of blue young stellar populations 
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is related to the strength of the SMBH accretion in our simulations. Therefore, Run 6 has 
the most significantly different color compared to other runs. 



Discussion and conclusions 



In this paper, we have investigated the effects of AGN feedback (in the form of com- 
bined radiative and mechanical energy and momentum deposition) on galaxies in different 
environments corresponding satellite and central galaxies. 

Although not including fully realistic stripping effects, our simulations predict various 
observable features that can be compared to the current and future observations, and which 
should change in a systematic way between central and satellite galaxies. The simulation re- 
sults need to be interpreted as a guide to understand a general consequence of environmental 
effects. Here, we discuss implications from our simulations. 

First of all, the local scaling relationship such as the mass ratio between the SMBH 
and its host galaxy can have a dispersion due to environmental effects. For example, if 
gas in a galaxy has experienced stripping, the growth of its SMBH is less disturbed than 
star formation which generally happens over the whole galaxy. Hence, the mass ratio be- 
tween the SMBH and the host galaxy of fixed galaxy mass could be significantly higher in 
satellite ellipticals in galaxy clusters. This also implies that a bias effect on deriving the 
scaling relationship probably exists due to selecting satellite galaxies less likely as observa- 
tion samples. Objects measured for their SMBH mass are generally bright objects which 
might not include galaxies having experienced gas stripping. Even though the r ecent mea- 



surements of the SMBH ni ass might have an intrinsic scatter of the ratio (e.g. iKim et al. 



2008bl : iGiiltekin et al.ll2009l ). the current measurements of the SMBH mass in satellite ellipti- 



cal galaxies have not been conducted systematically. Probably, the systematic investigation 



3e galaxies 
(Cote et al. 


m xne 
2006: 


Virgo Lviusxe 
Decarli et al. 


r migi 
2007; 


IX De useiui 
Gallo et al. 


CO xes 
2008) 



But we warn that the 

uncertainty and biases in measuring the SMBH mass might dominate the dispersions in the 



mass ratios even for galaxies in the Virgo cluster (e.g lBernardi et al.ll2007l ). 



The results of our simulations also imply the less frequent and earlier termination of 
AGN activity in satellite galaxies compared to central galaxies (i.e. galaxies with a large 
Rf). This result can be tested against the AGN fraction in galaxy clusters and groups (e.g. 
Gilmour et al.l 120071 ). The observed AGN fraction in galaxy clusters is generally higher than 
in field galaxies partially because of the high number density of galaxies (i.e. the effect from 
the cluster richness) and how AGN activities are selected in radio. X-ray, and IR observations 
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e.g. iHickox et al.ll2009l ). When considering this effect of the high nu mber density of galax- 



ies, A GN activity actually seems suppressed in galaxy clusters (e.g. iKoulouridis fc Plionis 
20101). The evolution of AGN fractions in galaxy clusters as a function of redshift also shows 
the increase in AGN activity at higher redshifts when satellite galaxies are still frequently 
active in their SMBH accretion, supporting our hyp othesis of environmental effects on AGN 
activity (IGalametz et al.ll2009l : iMartini et al.ll2009l ). The frequent multiple intensive SMBH 
accretion events in simulated central galaxies is consistent with several distinctive observable 
features of cluster central galaxies. For example, many central galaxies show multiple out- 
burst structures such as several gas clumps and surrounding ripples even though the origin 
of th ose observed struc t ures and what kind of AGN activity causes those is still un certain 
(e.g. 



Fabian et al.l 120051 : iGraham et al.l 120081 : iBlanton et al.ll2009t IClarke et al.l 120091 ) 



The very weak dependence of the ejected mass on Rf implies that satellite galaxies can 
be important sources of metal enrichment in the intracluster medium. Even though our 
simulation does not follow the evolution of the metal contents in the outflowed gas, the mass 
stripped out of the galaxy, which is about 2 x lO^^M© (see Figure H]), must be metal- rich being 
produced by gas recycling from the evolving stars and supernovae explosion. As explained in 
our models are similar to the leaky box model of galactic chemical evolution for satellites 



Pagell l2009l ) because of the outflowing boundary condition. Finally, the outflows from 
both satellites and central galaxies might cause mixi ng which is coupled to AGN activity 



and turbulence within the intracluster medium (e.g. lO'SuUivan et al.ll2005l ). The detailed 



abundance patterns in the cluster gas should thus reflect the outflow fr om the satellite 



cluster ellipticals in addition to the processes by a cluster central galaxy ( 


Domainko et al. 


2006; 


Schindler & Diaferio 


2008; 


Rasmussen & Ponman 


2009; 


Sivanandam et al. 


2009^. 



The recent star formation in cluster central galaxies (e.g. lRafferty et al.ll2008l : iKirkpatrick et al 



20091 ) is also naturally explained in our simulations by continuous recurrence of cooling flows 
and consequently heating by AGN feedback in galaxies with a large Rf. In our simula- 
tions, star formation from a cooling flow is temporally permitted while some amount of the 
cooling flow is accreted onto the SMBH and triggers AGN feedback which finally heats up 
and reverses the cooling flow as the cool gas is consumed by star formation. Contrary to a 
classical solution to the cooling fl ow problem in the cluster galaxies which stops the cooling 
flow completely (see lFabianlll994l . for a review), our solution to the cooling flow problem for 
the cluster central galaxies is simply permitting temporary rejuvenation of star formation 
and AGN activity which is self-regulated, leading to the termination of the cooling flow 



see 



Ciotti fc Ostrikerl 120071 . for a furt her discussiori) . This kind of solution can be called 



(Salome et al. 


2006; 


Bildfell et al. 


2008; 


Pipino et al. 



20091 : IWilman et al.ll2009l ). A further test of our results will be the Eddington ratio of the 
rejuvenated SMBH accretion in local cluster central galaxies because the SMBH accretion 
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rate in the simulated central galaxies is significantly lower than the Eddington accretion rate 
despite the change of Rf. 

To make our simulations more realistic, we need to improve our simulations by including 
two important physics: feedback by a radio jet from AGN activity and cosmological setup of 
environmental effects. In addition to the limits of our simulations already mentioned in our 
previous Paper II, for example, ID calculations, and possible different initial conditions such 
as models presented in Appendix, the radio jet is an important feedb ack component f or clus- 
ter c entral galaxies in environments of galaxy clusters and groups ( iBegelmanI |2004J : iKonigl 
20061 ). The observation of radio jets proves that the kinetic energy of jets is enough to be 
an effective feedback mode even though we do not have a clear explan ation about how much 
of th e relevant energy i s transported to ISM and intracluster medium (IMcNamara &: Nulsen 



20071 : IShin et al.l 1201 ll ). We note that powerful narrow jets such as that seen in the gi- 



ant elliptical M87 tend to drill through the ambient gas in the galaxy dep ositing some 



amount of energy with a relativistic fluid directly to the intergalactic medium (jFerrarilll998 



Owen et al.l 120001 ). Our implementation of stripping effect from different environments is 
simply parametrized by Rt in one dimension, and is not close to dynamically changing en- 



vironmental effects in cosm ological evolution (e.g. iTakeda et al.l Il984j : IVollmer et al.l 12001 



Toniazzo fc Schindlerll200l[ ). Moreover, our simulations do not include effects from tidal grav- 



itatio nal field which is commonly expected in cluster or group environments ( ID 'Er cole et al. 



2000|). Although the simple implementation makes the interpretation of simulation results 
obvious, the direct test of simulations will require better but complicated models of ram- 
pressure stripping and other environmental processes such as tidal stripping in galaxies which 
have realistic orbits in galaxy clusters and groups. 



A. Models with lower and higher velocity dispersions 

We also examine models with lower and higher initial stellar velocity dispersions of 
model galaxies than the fiducial model given in the main text. These models correspond 



(Ciotti et al. 


2009a: 


Shin et al. 


2010; 


Ciotti et al. 


2010) 



summarize results of the models Run 2, 4, 6, 8, and 10 with the lower velocity dispersion 
240 km/s and higher velocity dispersion 280 km/s than the fiducial velocity dispersion 260 
km/s. 

As shown in Figure [TT| the general trend of these models is not different from that of the 
fiducial model. One of the main differences is the frequency of the active SMBH accretion 
and star formation phase. As expected, models with the lower initial velocity dispersion 
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experience the less number of active phases than the others because of a less amount of 
stellar mass, which supplies materials for SMBH accretion and star formation. The second 
important feature is a high intensity of star formation rate in the beginning of the models 
with the higher velocity dispersion, while showing a lower SMBH accretion rate at the same 
time than the fiducial models. Due to a larger stellar mass in these models, the initial 
accumulation of cooling gas is large enough to trigger the intensive star formation in the 
beginning. But this intensive star formation competes with SMBH accretion, resulting in 
the lower SMBH accretion rate. In the late stage of evolution, the models with the higher 
velocity dispersion also show higher SMBH and star formation rates than the fiducial models. 
In particular, the SMBH accretion rate is close to the Eddington accretion rate. 

Figure [T2] shows that the general patterns found in Figure H] are still valid with the 
lower or higher initial stellar velocity dispersion. But the change of ratio between stellar 
and SMBH masses in the models with the higher velocity dispersion is less dependent of the 
different truncation radii. This weak sensitivity to the environmental effect is caused by the 
strong and more late SMBH accretion and the dominating high star formation rate in the 
models with the higher velocity dispersion. 

We are grateful to Michael Strauss, James Gunn, Gillian Knapp, Renyue Gen, Jenny 
Greene, and Ghristy Tremonti for useful discussions and careful reading. We thank the 
anonymous referee for comments which improved this manuscript. We also thank Yong Shi 
for the data used in Figure [2l M.-S.S. is supported by the Gharlotte Elizabeth Procter 
Fellowship of Princeton University. Gomputations were performed on the computational 
facilities of PIGSciE (Princeton Institute for Gomputational Science and Engineering). 
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Table 1: Truncation radius Rt of the hydrodynamical grid 



Run 


Radius (kpc) 


Run 


Radius (kpc) 


1 


51.4 


6 


239.7 


2 


88.5 


7 


287.3 


3 


127.2 


8 


344.4 


4 


166.9 


9 


377.0 


5 


218.9 


10 


412.8 



Note. — 



The effective radius of the galaxy models is 6.9 kpc at the beginning of the simulation. 
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O 




Fig. 1. — SMBH mass accretion rate and star formation rate. The SMBH mass accretion 
rate Mbh (blue) and star formation rate (red) show a strong time- dependence with mul- 
tiple intensive events consisting of several peaks. Mbh is usually lower than the Eddington 
accretion rate {short-long dashed line). Note how the frequency of the peaks of Mbh and 
depends on the truncation radius, particularly increasing from Run 1 {Rt = 51.4 kpc) to 
Run 3 {Rt = 127.2 kpc). 
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Fig. 2. — Star formation rate (M») and SMBH mass accretion rate (Mbh) around 3.23 Gyr 
{left top) and 13.31 Gyr {left bottom) and the distribution of M*/Mbh with respect to Mbh 
{right) in Run 6. Only early peaks of Mbh {thick solid line) are higher than the Eddington 
accretion rate {thick short-long dashed line) as shown in the left top panel. The early SMBH 
accretion is so strong that the following star formation is suppressed, while the late SMBH 
accretion is predominated by the concurrent intensive star formation. When M* {thin solid 
line) is high, Mbh is generally low in the early SMBH accretion as apparent in the left top 
panel. But when M* is high, Mbh is also high in the late SMBH accretion. The difference 
between the early SMBH accretion pattern and the late pattern is clearly found in the right 
panel showing the change for the time earlier than 3 Gyr {dashed line) and later than 7.5 
Gyr {solid line). Dots represent 3 and 7.5 Gyr, respec tively. Tr i angle s and triangles with 
arrows represent observed ratios and upper limits from lShi et al.l (120091 ). respectively. 
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Fig. 3. — Time evolution of mass components in the galaxy models. The different evolution 
of Mbh and for different values of Rt results in different final values of the mass ratios 
between the SMBH and stars, and different amounts of gas that are blown out. The mass 
AMbh added to the SMBH {thick solid line) is in general much smaller than the change 
in the stellar mass, AM^ {thin solid line). The total amount of X-ray emitting hot gas 
contained in the galaxy {thick short-long dashed line) decreases as Mbh, M^, and the total 
mass of outflowing gas {thin short-long dashed line) increase, particularly in the models with 
small truncation radii. 
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Fig. 4. — Comparison of the mass budget at 4, 8, and 14 Gyr for different values of Rt. The 
presentation follows the same line styles as given in Figure El At 2 Gyr differences among the 
simulations are not significant. As effects from different star formation and SMBH accretion 
history accumulate, the increase in stellar mass and the total mass of gas in Run 1 and 2 
begin to deviate from other simulations. But the amount of outflowing gas does not change 
much with respect to Rt- 



- 29 - 




Fig. 5. — Optical luminosity from accretion onto SMBHs and X-ray luminosity from hot 
ISM in their host galaxies for Run 1 and 6. In both simulations, the optical luminosity 
of SMBHs Lgy {thick solid line) correlates with the change of X-ray luminosity L^j^ {thin 
solid line). The optical luminosity from the accretion is usually below the 10% of the Ed- 
dington luminosity {short-long da shed line). We note tha t the optical luminosity is derived 
considering the extinction effects (ICiotti fc Ostrikerll2007l ). The bottom plots show the com- 
plex evolution of Run 1 at about 3.9 Gyr and of Run 6 at about 13.3 Gyr with bolometric 
luminosity from SMBH without dust extinction {dotted line). When L'^^, i.e. the SMBH 
accretion rate, is high, shows the oscillatory structure because of feedback from the 

accreting SMBH. Yet, the accumulated cold gas finally strongly extincts optical luminosity 
from the accreting SMBH while -Z^iIm SMBH bolometric luminosity increase around 

3.86, 3.92, and 13.31 Gyr. 
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Fig. 6. — Radial distribution of the ISM temperature and density in Run 1 and 6. The 
truncation radii are indicated as dashed hues. At 14 Gyr {thin line), Run 1 and 6 are found 
in the low-luminosity hot accretion phase; the higher central temperature and density in Run 
6 are because of the influence of the larger final SMBH mass. However, Run 6 still shows 
the effect from its recent active phase forming the sharp shell structure. When both models 
are in the active phase at 3 Gyr [thick line) either increase of core temperature or shell 
structures appears in the radial distributions. Triangles, squares, and circles correspond to 
observed values at 2 an d 10 kpc in NGC 4125, NGC 720, and NGC 6482, respectively, from 



Humphrey et al.l (120061 ). 
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Fig. 7. — Episodic and net lifetime of optically luminous quasar phase (i.e. Mb < —23 
mag) in Run 1, 4, 6, and 9. The episodic hfetime is presented for each peak of Mbh as cross 
symbols. The net lifetime is given in the parenthesis next to the simulation number. The 
episodic lifetime increases as the simulation continues until about 3.5 Gyr when the peaks 
of SMBHs accretion rate are higher than the Eddington accretion rate. 
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Fig. 8. — Outflow velocity and density as a function of Rt at 4, 8, and 14 Gyr. Velocity 
{solid line) and density {dashed line) is estimated on the last radial grid point Rf which 
represents the boundary of gas truncation. Models with small Rt have the more dense and 
fast flows that pass out the truncation radius, producing the amount of mass loss which is 
comparable to that of models with the large Rt when we measure the mass loss at 10ii!e- We 
note that the mass flowing out at late times is neghgible for large values of Rt- 
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Fig. 9. — Systematic color-magnitude diagram at 8.5 Gyr {dashed line) and 14 Gyr [solid 
line) without the optical radiation from the SMBH. From the star formation history of eac h 
run, we retrieve the synthesized spectra by using the BC03 model (IBruzual &: Charlotll2003l ). 
Solar metallicity is assumed for black lines, while half solar metallicity is used for blue lines. 
The effects from dust extinction are not considered in the construction of spectra. The 
numbers in circles represent the name of the run. The arrows represents color excess for 
Ar = 0.1 with the dust e xtinc t ion cu rves of the Milky Way [solid line) and Small Magellanic 
Cloud [dashed line) from iPeil (jl992[ ). respectively. 
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Fig. 10. — Same as Fig. [9l but with the optical radiation from the SMBH. The include d 
quasar spectrum is a luminosity-scaled spectrum of type-I quasar (IVanden Berk et al.ll200ll ). 
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Fig. 11. — Same as Fig. [H but with the initial velocity dispersions of 240 km/s (left) and 
280 km/s (nghi) for Runs 2, 4, 6, 8, and 10. 
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Fig. 12. — Same as Fig. HI but with the initial velocity dispersions of 240 km/s (left) and 
280 km/s (nght) for Runs 2, 4, 6, 8, and 10. 



